Deflation of Eigenvalues for Iterative Methods in Lattice QCD
نویسندگان
چکیده
Work on generalizing the deflated, restarted GMRES algorithm, useful in lattice studies using stochastic noise methods, is reported. We first show how the multi-mass extension of deflated GMRES can be implemented. We then give a deflated GMRES method that can be used on multiple right-hand sides of Ax = b in an efficient manner. We also discuss and give numerical results on the possibilty of combining deflated GMRES for the first right hand side with a deflated BiCGStab algorithm for the subsequent right hand sides.
منابع مشابه
Deflation of Eigenvalues for GMRES in Lattice QCD
This paper looks at the iterative solution of complex, non-Hermitian systems of linear equations associated with a lattice QCD problem in particle physics. Let the n by n system of equations be Ax = b. We use a version of the GMRES method that deflates eigenvalues and can solve systems with multiple right-hand sides. We give examples that show deflating eigenvalues can make a significant improv...
متن کاملLocal coherence and deflation of the low quark modes in lattice QCD
The spontaneous breaking of chiral symmetry in QCD is known to be linked to a non-zero density of eigenvalues of the massless Dirac operator near the origin. Numerical studies of two-flavour QCD now suggest that the low quark modes are locally coherent to a certain extent. As a consequence, the modes can be simultaneously deflated, using local projectors, with a total computational effort propo...
متن کاملRestarted Block Gmres with Deflation of Eigenvalues
Block-GMRES is an iterative method for solving nonsymmetric systems of linear equations with multiple right-hand sides. Restarting may be needed, due to orthogonalization expense or limited storage. We discuss how restarting affects convergence and the role small eigenvalues play. Then a version of restarted block-GMRES that deflates eigenvalues is presented. It is demonstrated that deflation c...
متن کاملComparing iterative methods to compute the overlap Dirac operator at nonzero chemical potential
The overlap Dirac operator at nonzero quark chemical potential involves the computation of the sign function of a non-Hermitian matrix. In this talk we present iterative Krylov subspace approximations, with deflation of critical eigenvalues, which we developed to compute the operator on large lattices. We compare the accuracy and efficiency of two alternative approximations based on the Arnoldi...
متن کاملDeflation acceleration of lattice QCD simulations
Close to the chiral limit, many calculations in numerical lattice QCD can potentially be accelerated using low-mode deflation techniques. In this paper it is shown that the recently introduced domain-decomposed deflation subspaces can be propagated along the field trajectories generated by the Hybrid Monte Carlo (HMC) algorithm with a modest effort. The quark forces that drive the simulation ma...
متن کامل